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Abstract 
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Comparing perturbative calculations with a lattice computation of the static energy in quantum chromodynamics at 
short distances, we obtain a determination of the strong coupling a s . Our determination is performed at a scale of 
around 1.5 GeV (the typical distance scale of the lattice data) and, when evolved to the Z-boson mass scale Mz, it 
corresponds to a s (M z ) = 0.11561° ™^. 
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The energy between a static quark and a static an- 
tiquark that are separated a distance r, i.e. the quan- 
tum chromodynamics (QCD) static energy, Eo(r), is a 
good object to study in order to understand the behav- 
ior of the theory. One can identify a long-distance part 
and a short-distance part of the static energy, both of 
which can be computed with lattice simulations. Here 
we will focus only on the short-distance part, i.e. on 
distances r < 0.234 fm, where weak-coupling calcu- 
lations are also reliable. The comparison of the lattice 
computation with the perturbative calculation tests our 
ability to describe the short-distance regime of QCD, 
and provides information on the region of validity of 
the weak-coupling approach. This comparison also al- 
low us to determine the strong coupling a s , which is the 
subject of this talk. There has been a lot of recent ac- 
tivity regarding both, the lattice computations and the 
perturbative calculations of the static energy, which has 
allowed for a precise determination of or, from it to be 
possible 111]. 

On the perturbative side, the static energy is known, 
at present, including terms up to order aq +n In" a s with 
n > o M. That is, three-loop with resummation at 
sub-leading accuracy of the In as terms that appear at 
short distances. We refer to this level of accuracy as 
next-to-next-to-next-to leading-logarithmic (N 3 LL). 

On the lattice side, the static energy has recently been 
calculated in 2 + 1 flavor QCD H, using a combina- 
tion of tree-level improved gauge action and highly- 
improved staggered quark action [10]. This computa- 
tion employed the physical value for the strange-quark 
mass m, and light quark masses equal to m s /2Q, which 



correspond to a pion mass of about 160 MeV in the con- 
tinuum limit, very close to the physical value. The com- 
putation was performed for a wide range of gauge cou- 
plings, and was corrected for lattice artifacts. It allows 
to study the static energy down to distances r 0.065 
fm. 

The perturbative expressions for the static energy de- 
pend on the value of the QCD scale A^g- (in the MS 
scheme), and we can use the comparison with lattice 
data to determine it. For that, we assume that pertur- 
bation theory (after implementing a cancellation of the 
leading renormalon singularity) is enough to describe 
lattice data in the range of distances we are consider- 
ing. Then, the general idea is that we can search for the 
values of A^g- for which the agreement with lattice im- 
proves when the perturbative order of the calculation is 
increased; and in that way find the values of A^ that 
are allowed by lattice data. This same program was al- 
ready performed for the quenched case in Ref. the 
unquenched computation of the static energy in Ref. |0] 
allows us to do the same here in the unquenched case, 
and therefore obtain a value for a s . 

The static energy on the lattice is calculated in units 



of the scales ro or r\, defined as 111 1111211 
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we use the values of ro or r\ in Ref. Igy to obtain A^g- in 
physical units. In the perturbative calculation one needs 
to implement a scheme that cancels the leading renor- 
malon singularity This kind of schemes introduce 



1 In the lattice computation the results calculated at different lattice 
spacings are normalized to a common value at a certain distance. 
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an additional dimensional scale in the problem (that we 
denote as p). We implement the renormalon cancel- 
lation according to the scheme described in Ref. 11411 : 
then, the natural value of the scale p is at the center of 
the range for which we have lattice data. But since any 
value of p around this natural value cancels the renor- 
malon, we can exploit this freedom to search for a set 
of p values that allow for an optimal determination of 
roAj^g. To obtain our central value for roA^o we ^ et 
p vary around its natural value; then, for each value of 
p and at each order in the perturbative expansion, we 
perform a fit to the lattice data (roA^g- is the parame- 
ter of the fits); and finally select the p values for which 
the reduced y 1 of the fit decreases when increasing the 
perturbative order. Our central value for r^K-^ is then 
given by the average (weighted by the inverse x 2 ) of 
those fit valued We can perform the above analysis at 
different orders of accuracy; at N 3 LL accuracy the per- 
turbative expression depends on an additional constant 
(due to the structure of the renormalization group equa- 
tions at this order), which would also need to be fitted to 
the lattice data. When we do the fits at N 3 LL accuracy, 
we find that the x 1 as a function of r^K-^ is very flat; 
we interpret this fact as the data not being sensitive to 
the sub-leading ultrasoft logarithms. Therefore we take 
the result at three loops with resummation of the leading 
ultrasoft logarithms as our best result. For illustration, 
we present here the expressions for the static energy at 
this level of accuracy. The static energy at next-to-next- 
to-next-to leading order (N 3 LO) is given by 
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2 Note that the absolute value of the x 2 (from the fits of the pertur- 
bative expressions to the lattice data), at a given order in perturbation 
theory, does not have any particular significance. The reason for that 
is that, in principle, one does not know exactly how large the terms at 
the next perturbative order will be. Therefore, one does not know how 
accurately the perturbative expression should describe the lattice data 
at a given order. For that reason, our procedure to determine roA^-, 
described above, does not use absolute values of the ^ 2 s, but rather 
comparisons between x 2 values at different perturbative orders. Oth- 
erwise there would be the danger that one artificially reduces the x 2 
at a low perturbative order, by using a "wrong" value of roA^r that 
is not suitable for a more precise expression at a higher perturbative 
order. 



with (numerically, for N c - 3) 

a, = 23.032 - 1.880771/, (3) 
a 2 = 1396.3- 192.9071/ + 4.9993n}, (4) 
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-20.6009n}, 



(5) 



where C F = (N^ - l)/(2N c ), C A = N c , n f is the number 
of light flavors (i.e. «/ = 3 in our case), and K\ is a con- 
stant that, in the comparison with data, gets absorbed 
in the constant used to make the static energy coincide 
with the lattice point at the shortest distance available. 
If we include the resummation of the leading ultrasoft 
logarithms we have 
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where p is the ultrasoft scale (of order a s /r), and fio = 
(U/3)C A - (4/3)7Vn/, with 7> = 1/2. To imple- 
ment the required renormalon cancellation we use the 
so-called RS scheme II 1411 . That is, if we calculate the 
static energy at m-loop order in perturbation theory we 
add the following term to it Ja] 

RSsubtr. = R'PYXte) as(pT+1 

n=\ 

Z T (n-i-l+b-1c) 
J 1 ni+»-» ■ <7) 

with R s = -1.123 the normalization of the u = 1/2 
renormalon singularity (which we computed approxi- 
mately according to the procedure in Ref. lfl5ll ). and 

do - 1 , 
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(the higher order coefficients of the beta function, ^61,2,3, 
can be found, for instance, in Refs. |T^, 171). 

Having determined our central value for roArr? we 
now need to assign an error to it. The error must re- 
flect the uncertainties associated to the neglected higher- 
order terms in the perturbative expansion. To account 
for that, we consider the weighted standard deviation 
in the set of p values we found before, and the differ- 
ence with the weighted average computed at the previ- 
ous perturbative order. The latter term turns out to be 
the dominant error; we then add the two errors linearly. 
Additionally, we also redo the analysis with alternative 
weight assignments (p-value, and constant weights); we 
obtain compatible results, and quote and error that cov- 
ers the whole range spanned by the three analyses. As 
a further cross-check, we can compare the analysis per- 
formed with the static energy normalized in units of ro 
(our default choice) and the one with the static energy 
normalized in units of r\ ; we find that the two analyses 
give consistent results^. 

Our final result reads 



r Q K m = 0.70 + 0.07, 



(10) 



which using the value of ro from Ref. ||9|] gives 

a s (p= 1.5GeV,n / = 3) = 0.326 ±0.019. (11) 

When we evolve Eq. (fTTT i to the Z-mass scale, M z , we 
obtain 



a s (M z ,n f =5) = 0.1156to 
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where we have used the Mathematicapackage RunDec 
111 811 to obtain the above number (4 loop running, with 
the charm quark mass equal to 1 .6 GeV and the bottom 
quark mass equal to 4.7 GeV). 

We compare our result with other recent lattice deter- 
minations of el's in Fig.Q] Our central value is a bit lower 
than those of the other lattice determinations. 

Our determination is performed at a scale of around 
1.5 GeV. This scale corresponds to (the inverse of) the 
typical distance where: (i) we have lattice data, and (ii) 
the weak-coupling calculation is reliable. This means 
that our analysis represents the lowest-energy determi- 
nation of or s available, and can therefore be an important 
new ingredient to further test the running of a s . Previ- 
ously, the lowest-energy determination was that com- 
ing from hadronic t decays (performed at m T = 1.78 
GeV). For comparison, the value of the pre-average of 



'This is a cross-check and not just a trivial reseating, because the 
errors and systematics entering the lattice analyses normalized in units 
of ro or n are different. 
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Figure 1 : Comparison of our result (red -lighter- point) with other re- 
cent lattice determinations of a s (blue -darker- points). The references 
are: HPQCD f 19], JLQCD 01, PACS-CS ETM 01 



a s determinations from r decays that is currently used 
by the Particle Data Group (PDG) Q is a s (M z ) = 
0.1197 + 0.0016. 

To summarize, we have obtained a determination of 
or s by comparing perturbative calculations with a lat- 
tice computation of the short-distance part of the QCD 
static energy. Our determination is at three-loop ac- 
curacy (including resummation of the leading ultrasoft 
logarithms), and is performed at a scale of 1 .5 GeV (and 
therefore it constitutes the lowest-energy determination 
of a s available). When evolved to the scale M z , it cor- 
responds to a s (M z ) — 0.1 156^qqq22- A very important 
outcome of our work is also that our analysis shows, for 
the first time in QCD with nj = 2 + 1 flavors, that per- 
turbation theory (after cancellation of the leading renor- 
malon singularity) can describe the short-distance part 
of the static energy. This is illustrated in Fig. [2] 
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